Net-proton probability distribution in heavy ion collisions 
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We compute net-proton probability distributions in heavy ion collisions within the hadron reso- 
nance gas model. The model results are compared with data taken by the STAR Collaboration in 
Au-Au collisions at y/SNN ~ 200 GeV for different centralities. We show that in peripheral Au-Au 
collisions the measured distributions, and the resulting first four moments of net-proton fluctuations, 
are consistent with results obtained from the hadron resonance gas model. However, data taken 
in central Au-Au collisions differ from the predictions of the model. The observed deviations can 
not be attributed to uncertainties in model parameters. We discuss possible interpretations of the 
observed deviations. 
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I. INTRODUCTION 

One of the objectives of heavy ion experiments at 
CERN and BNL is to probe properties of the QCD phase 
diagram related to deconfinement and chiral symmetry 
restoration. To experimentahy verify the phase change 
in a medium created in such collisions, one needs observ- 
ables that are sensitive to critical behavior. Modifications 
in the magnitude of fluctuations or in the corresponding 
susceptibilities of conserved charges have been suggested 
as possible signals for chiral symmetry restoration and 
deconfinement ^^■ 

Fluctuations of baryon number and electric charge 
diverge at the hypothetical critical end point in the 
QCD phase diagram at non-zero temperature and baryon 
chemical potential, while they remain finite along the 
cross-over boundary. Consequently, large fiuctuations of 
baryon number and electric charge as well as a non mono- 
tonic behavior of these fiuctuations as function of the col- 
lision energy in heavy ion collisions have beenproposed 
as a signature for the QCD critical end point [E H 0] • 

It also has been argued that even in the absence of a 
critical end point, fiuctuations of conserved charges, and 
their higher order cumulants, can be used to identify the 
phase boundary. It is expected that the fiuctuations are 
modified if the chemical freeze-out, i.e., the generation 
of hadrons and their fiuctuations, occurs shortly after 
the system passed through a region where quarks were 
deconfined and chiral symmetry was partially restored 
[a, HjB- Thus, fiuctuations in the hadronic phase can 
show structures attributed to critical dynamics at the 
restoration of chiral symmetry. 

The critical region characterizing the cross-over transi- 
tion in the QCD phase diagram is expected to be located 
close to the freeze-out curve extracted from heavy ion 
experiments [9| . On this phenomenologically determined 



curve all particle yields achieve their measured values 
jlO|-|l^. Thermodynamics at freeze-out is, to a first ap- 
proximation, well described by the hadron resonance gas 
(HRG) model, which was shown to be very successful in 
describing not only the data on particle yields but also 
the thermodynamics of a strongly interacting medium at 
low temperature as computed in lattice QCD |3. Il4| - KL8| . 

If chemical freeze-out occurs near or at the QCD phase 
boundary, this should be reflected in the higher order 
cumulants of charge fluctuations since the sensitivity to 
critical dynamics grows with increasing order [7[ . Conse- 
quently, the values of higher order cumulants can differ 
significantly from the results of the HRG along the freeze- 
out curve even if lower order cumulants agree. In particu- 
lar, at vanishing chemical potential, the sixth and higher 
order cumulants can even be negative in the hadronic 
hase while the HRG yields positive values everywhere 
In fact, the results of the HRG model on cumulants 
of charge fluctuations can serve as a theoretical baseline 
for the analysis of heavy ion collisions [2|, 0, [ij] . In equi- 
librium, any deviation from the HRG model would be a 
reflection of genuine QCD properties not accounted for 
by the model and could constitute evidence for critical 
phenomena at the time of hadronization. 

Recently first data on charge fluctuations and higher 
order cumulants, identifled through net-proton fluctua- 
tions, were obtained by the STAR Collaboration in Au- 
Au collisions at several collision energies [1^ [20] . To ex- 
plore possible signs of criticality, the STAR data on the 
flrst four cumulants were compared to HRG 0, 120] and 
lattice QCD [2l|, [22[ results. The basic properties of the 
measured fluctuations and ratios of cumulants are con- 
sistent with the expectations based on HRG as well as on 
lattice QCD calculations. However, a more detailed com- 
parison of the HRG model with STAR data reveals that 
at high energies, deviations can not be excluded [7l.l2C|. 
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FIG. 1: Left-hand figure: Net-proton distributions calculated in the hadron resonance gas using Eq. |4] and compared with 
STAR data for Au-Au collisions at \/sjviv = 200 GeV for different centralities [ly]. The full- lines are obtained with experimental 
inputs for proton (Np) and anti-protons {Np) yields, while the broken-lines are obtained with (Np) and (Np) computed in the 
thermal model with the parameters at chemical freeze-out taken from Ref. [13|. Right-hand figure: Mean (M), variance (a), 
skewness (S) and kurtosis (k) calculated from the probability distributions shown in Fig. [T] (left) for different centralities. Data 
are from STAR fTl. 



All moments of net-proton fluctuations as well as the 
related cumulants can be calculated once the underlying 
probability distribution is known. Therefore, it is inter- 
esting to confront the net-proton distributions obtained 
for the HRG with those measured in heavy ion collisions. 
Since the probability distributions contain information 
on all cumulants, such a comparison may provide useful 
insights into the origin of possible deviations from the 
HRG baseline. An analysis of this kind may also provide 
additional information on the relation between chemical 
freeze-out and the QCD cross-over transition or even on 
the existence of a critical end point in the QCD phase 
diagram. 

In the following we calculate the net-proton probability 
distributions in the HRG model. We show that in this 
model the net-proton distribution can be expressed solely 
in terms of the measurable yields of protons and anti- 
protons. The HRG model results are compared with data 
taken by the STAR Collaboration in Au-Au collisions at 
RHIC at y'sNN = 200 GeV. For central collisions at this 
energy, the HRG model yields a distribution which is 
broader than the one seen experimentally. We discuss the 
possible origin of the deviations and point out that they 
are in qualitative agreement with the differences between 
HRG model calculations and (lattice) QCD results found 
in the vicinity of the QCD phase boundary and above. 



II. PROBABILITY DISTRIBUTION OF 
CONSERVED CHARGES 



Consider a sub-volume y of a thermodynamic system 
described by the grand canonical ensemble consisting of 
charged particles q and anti-particles g at a given temper- 
ature T and chemical potential /i. The latter is related 
to the conserved net charge N — Nq — Nq. The prob- 



ability distribution P{N) for finding a net charge num- 
ber N in the volume V is given in terms of the canoni- 
cal Z(T, V, N) a nd g rand canonical Z(T, T/,/i) partition 
functions [iSllllll, 



P{N) = Z{T,V,N)e^''^-^^P^^'^'^ . 



(1) 



Here we have introduced the thermodynamic pressure, 
InZ = Vl3p{T, ^), as well as the shorthand notation jl = 
f3fi and /3 = 1/T. 

The canonical partition function Z{T, V, N) can be ob- 
tained from the thermodynamic pressure through the dis- 
crete Fourier transform 



Z{T,V,N) 



1 

2^ 



-t^NBVp{T,i^/f3) 



(2) 



where the chemical potential was Wick rotated by the 
substitution /t — > i(j). Eqs. ([1]) and © define the proba- 
bility distribution of a conserved charge in a sub- volume 
I^ of a thermal system described by the thermodynamic 
pressure p(T, ^). 

In the following we focus on fluctuations of the net- 
baryon number and consider the corresponding proba- 
bility distribution in a strongly interacting medium. We 
model the thermodynamics of this system using the HRG 
partition function, which contains all relevant degrees of 
freedom in the hadronic phase and implicitly includes 
the interactions responsible for resonance formation. In 
this model the thermodynamic pressure is a sum of me- 
son and baryon contributions. This implies that only the 
baryonic pressure pb contributes to the probability dis- 
tribution of the net-baryon number. In the HRG model 
Pb consists of contributions from all baryons and bary- 
onic resonances (lO[. In the Boltzmann approximation, 
(3VpBiT, n) = b{T, V, n) -h b{T, V, n), where b and b is the 



mean number of baryons and anti-baryons respectively, 



27r2 

z^baryons 

Here qi = {Si,Qi) is a two-component vector composed 
of the strangeness and electric charge carried by particle 
i, flq = {^Ji.s,IJ^Q) the corresponding chemical potential 
vector, gi the spin-isospin degeneracy factor and K2 is 
a modified Bessel function. The mean number of anti- 
baryons b is obtained by the substitution /z — >■ — /z for all 
relevant chemical potentials. 

In the hadron resonance gas model, the canonical 
partition function Z{T,V, N) is computed directly from 
Eq. ([2Jj_^using the thermodynamic pressure discussed 
above [l3] ■ The resulting probability distribution of net- 
baryon number can be expressed solely in terms of the 
mean number of baryons and anti-baryons, 



P{N) 



N/2 



In{^> 



cxp[-{b + b)] 



(4) 



where /Ar(x) is a modified Bessel function. 

We note that the above arguments remain valid also for 
subsystems which are limited not only in position space, 
but more generally in phase space. In particular, the 
introduction of cuts in momentum space leave Eqs. ([l}, 
([2| and ([4]) unchanged after an appropriate redefinition of 
the partition functions and densities. Moreover, the re- 
striction to one particle species, e.g. protons with proper 
account for resonance decays, is easily accommodated. 



III. PROBABILITY DISTRIBUTION OF THE 
NET-PROTON NUMBER 

In this section we confront the HRG model results 
for the probability distribution of the net-proton num- 
ber with data of the STAR Collaboration [l^, ^. The 
data is obtained at mid-rapidity in a restricted range of 
transverse momentum, 0.4 GeV < pr ^ 0.8 GeV. The 
probability distribution for protons is readily obtained 
from Eq. ^, by replacing the mean number of baryons 
b and anti-baryons b by that of protons (Np) and anti- 
protons (Np), respectively. 

Using Eq. 1^ we can compute the net-proton distri- 
bution provided we have access to the mean values {Np) 
and (Np) measured in the same kinematic window. For 
Au-Au collisions at ^/snn = 200 GeV, STAR data on 
the pf-distribution of protons and anti-protons is avail- 
able at several centralities [l3|. By integrating the pt- 
spectra of anti-protons in the pt-window where the net- 
proton number was obtained, we find: (Np) = 5.233(95), 
(Np) = 1.838(7) and {Np) = 2.844(98) for (0-5%)-central, 
(30-40)%-mid-central and (70-80)%-peripheral colhsions, 
respectively. Since the proton data were measured in 
a slightly larger pt window, we avoid systematic errors 
that may arise from an extrapolation by computing (Np) 



from the net-proton yields {M — (Np) — {Np)), with 
M ~ 1.715, M ~ 0.597 and M ~ 0.08 for central, mid- 
central and peripheral collisions [1^ . 

In Fig. [T] (left) we compare the STAR data on the net- 
proton multiplicity distribution in Au-f-Au collisions at 
y^sj^ =200 GeV [H [13 with the probability distribu- 
tions obtained in the HRG model (Eq. (j4])), using the 
experimental data on M and (Np) as input. The data 
correspond to several centrality bins in the rapidity win- 
dow \y\ < 0.5. The distribution in Eq. ^ is normalized 
to unity. In order to confront the HRG model with data 
on an absolute scale, we adjust the normalization to that 
of the experimental data in each centrality bin. 

Fig. [1] (left) shows that in peripheral collisions at 
^/sNN =200 GeV the HRG model results agree well with 
the data. However, with increasing centrality deviations 
develop; in central collisions the hadron resonance gas 
yields a distribution, which is broader than the exper- 
imental one. Below we argue that such deviations are 
expected if the freeze-out conditions probed by fluctu- 
ations are located close to the QCD cross-over temper- 
ature. Consequently, the deviations observed in Fig. [T] 
(left) and the dependence on centrality, could be an in- 
dication for critical behavior. 

We have also calculated the number of protons and 
anti-protons with chemical freeze-out parameters ad- 
justed to the multiplicities measured in Au-Au collisions 
at y/SNN =200 GeV. The centrality dependence of the 
freeze-out parameters were determined in Ref. [l3|. In 
this approach an additional input parameter is required, 
the effective volume (V), which we fix by requiring that 
the measured net-proton number M at chemical freeze- 
out is reproduced by the HRG model. As shown in 
Fig. [1] (left) the probability distributions obtained in this 
way are consistent with those computed directly from 
the measured yields of protons and anti-protons using 
Eq. Qj. The consistency of the two approaches strength- 
ens our conclusion that the HRG model does not describe 
the net-proton probability distribution P{N) for central 
Au-Au collision at ^/sNN ~ 200 GeV. 

Recently the STAR Collaboration also presented pre- 
liminary data on the net-proton distribution in Au-Au 
collisions at ^/snn = 39 GeV for various centralities [201 ■ 
It would be interesting to compare this data with the 
HRG model. In this case, however, the corresponding 
data on pi-distributions, which would allow one to deter- 
mine (Np) and (Np) directly from the experiment is not 
available. Thus, here we can only follow the second ap- 
proach, i.e., employ the measured net-proton number M 
and the chemical freeze-out parameters from Ref. |12| to 
determine (Np) and (Np), which are then used as input 
for the calculation of P{N) based on Eq. (jj]). 

A first comparison of the HRG model calculation 
with the probability distribution obtained for central 
Au-Au collisions at ^/sNN = 39 GeV indicates that at 
this energy the shape and magnitude of the measured 
net-proton distribution, and consequently also the first 
four measured moments, are described well by the HRG 



model. This suggests that in central Au-Au collisions at 
y^sNN = 39 GeV the fluctuations as well as the particle 
yields are characterized by the thermodynamic freeze-out 
conditions corresponding to the statistical operator of the 
hadron resonance gas model. Clearly, this result needs 
to be confirmed by the final STAR data at this lower en- 
ergy. If correct, this would suggest that the deviations 
from the HRG results grow with increasing energy. This 
feature can also be tested at the LHC. 

The observed deviations in the probability distribution 
are also manifested in differences between calculated and 
measured cumulants of the net-proton fluctuations. In 
Fig. [1] (right) we show the mean, variance, skewness and 
kurtosis obtained from the probability distributions in 
Fig. [T] (left). The HRG model, with experimental input 
for proton and anti-proton yields a slightly better de- 
scription of all four moments. However, a good overall 
description is obtained only for peripheral collisions. 

The fact, that the HRG model yields a distribution 
which is broader than the experimental one implies that 
deviations arise already on the level of the second order 
cumulant (variance), which has the smallest experimen- 
tal error. This is expected if the particle freeze-out occurs 
near the QCD cross-over transition. In the cross-over re- 
gion the baryon number susceptibility (xf ), i-e., the sec- 
ond order cumulant a'^ = VT^X2 j keeps rising steeply 
with temperature in the HRG model, while in QCD cal- 
culations it bends over and eventually approaches a finite 
value at high temperatures. In the Gaussian approxi- 
mation (the leading order cumulant expansion) to the 
probability distribution, P{N) ^ cxp [-7X^^(20"^)], this 
implies that the distribution in QCD is narrower than in 
the HRG model. 

A possible interpretation of this effect may be related 
to the proximity of the freeze-out and cross-over regions 
probed at the highest beam energy. Lattice calculations 



suggest that at larger /Lt, the freeze-out curve the cross- 
over transition separate [J, [l^. Hence, one expects that 
the fluctuations reflect the critical dynamics at the cross- 
over transition only at the highest beam energies. 



IV. CONCLUSIONS 

We have analyzed properties of the net-proton prob- 
ability distributions in heavy ion collisions within the 
hadron resonance gas model. In this model these dis- 
tributions can be expressed solely in terms of the mean 
numbers of protons and anti-protons in a thermal system. 
This provides a direct and unambiguous way to compare 
experimental data with model predictions. 

We have shown that the HRG model describes 
the data obtained by STAR in Au-Au collisions at 
y^SNN — 200 GeV for peripheral events but differ for 
central ones. Since the probability distributions are com- 
puted directly from the measured proton and anti-proton 
yields, this deviation is not due to uncertainties in model 
parameters. We suggest that this effect could be due to 
the proximity of the freeze-out and cross-over regions at 
the highest beam energy. In order to substantiate this 
interpretation, data on the net-proton distribution at 
lower RHIC energies and at LHC energies are needed. 
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